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PREFACE 


This  report  describes  the  modifications  made  to  the  Air  Force  MODTRAN2  program  in 
order  to  create  the  software  model  AURIC-M.  Input  instructions  for  AURIC-M  are  also  provided, 
as  a  supplement  to  the  MODTRAN  User's  Manual.  AURIC-M  resulted  from  close  technical 
collaboration  between  Aerodyne  Research,  Inc.  and  the  Phillips  Laboratory  Geophysics 
Directorate,  in  particular  L.A.  Hall  (the  contract  monitor),  Gail  Anderson  and  Jim  Chetwynd. 
Special  acknowledgement  is  made  of  the  contributions  of  Ken  Minschwaner  (New  Mexico  Tech.), 
who  developed  the  Schumann-Runge  model  and  who  consulted  to  Aerodyne  Research  on  the 
program.  Contributors  at  Aerodyne  Research  were  Curt  Betchly  (now  at  Lockheed, 

Ft.  Worth  Corp.),  Robert  Huguenin  (now  at  Applied  Analysis,  Inc.),  Michael  Weinberg, 

Malcolm  LeCompte,  and  Robert  Hickey. 
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1.0  INTRODUCTION 


The  ultraviolet  (UV)  spectral  region  has  current  strong  interest  for  both  military  and  pure 
scientific  studies.  Military  research,  for  example,  often  utilizes  the  UV  emissions  from  rocket 
motor  exhausts,  and  the  possibility  for  smaller  optical  sensors  than  for  the  infrared  (IR).^ 
Wavelengths  below  300  nm  are  important  because  atmospheric  ozone  and  oxygen  strongly 
attenuate  the  solar  irradiance,  leaving  a  dark  background  for  target  detection. 

A  major  area  of  scientific  interest  is  in  atmospheric  photochemistry,  in  which  the  UV  region 
is  of  prime  importance.23  The  absorption  of  O2  and  O3  largely  determine  the  extent  of  penetration 
of  sunlight  into  the  stratosphere  and  troposphere,  where  it  causes  photodissociation  of  a  variety  of 
molecules.  The  photodissociation  of  O2  is  especially  important  as  the  major  source  of  odd 
oxygen.  3 

The  Air  Force  Phillips  Laboratory,  Geophysics  Directorate,  has  long  had  a  successful 
program  for  the  development  of  computer  models  of  atmospheric  transmission  and  radiance, 
including  the  widely-used  LOWTRAN  and  MODTRAN  series.^-^  These  models  have  strong 
capability  in  the  infrared,  and  moderate  capability  in  the  visible  and  UV.  The  ongoing  development 
of  AURIC,  the  Atmospheric  Ultraviolet  Radiance  Integrated  Code,  will  provide  a  similar  capability 
with  specific  attention  to  the  UV.l 

One  aspect  of  the  AURIC  development  is  the  improvement  of  MODTRAN2  in  the  UV, 
resulting  in  a  model  called  AURIC-M.  This  Fortran  program  extends  MODTRAN2  both  spectrally 
(out  to  83,(KX)  cm*^  or  120  nm)  and  in  altitude  (up  to  1000  km),  and  provides  improved  solar 
irradiance  and  O2  absorption  models. 

In  addition,  AURIC-M  can  compute  molecular  photodissociation  rates,  including  those  of 
chlorofluorocarbons  (CFCs).  This  report  describes  the  development  of  AURIC-M,  and  provides 
instructions  for  its  use.  Section  2  contains  the  algorithmic  and  physical  model  descriptions,  while 
Section  3  describes  the  changes  in  software  usage  from  that  of  MODTRAN.  The  references  are 
listed  in  Section  4. 
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2.0  TECHNICAL  DEVELOPMENT 


The  AURIC-M  model  was  developed  from  the  existing  MODTRAN  model.^  Changes 
were  made  in  a  number  of  areas,  including: 

1 .  Inclusion  of  better  O2  line  and  continuum  absorption  models  from  40,600  cm-1  to 

66,000  cm-1, 

a)  Replacement  of  the  older  Schumann-Range  band  model  with  a  temperature- 
dependent  model  from  49,000  cm-l  to  57,000  cm-^  at  1  cm'^  intervals, 

b)  Replacement  of  the  older  Herzberg  continuum  model  above  40,600  cm-^ 
with  a  new  Herzberg  and  Schumann-Range  combined  continuum  to 
66,000  cm-^ 

2 .  Increase  in  the  number  of  atmospheric  layers  available  for  doing  the  calculations, 
easily  changed  with  a  re-compilation, 

3 .  Increase  in  the  altitude  extent  of  the  atmospheric  profile  data  for  Model  6,  the 
1976  US  Standard  Atmosphere,  to  1000  km, 

4.  Increase  in  the  spectral  extent  and  resolution  of  the  solar  irradiance  model  using 
1987  SUSIM  data, 

5 .  Inclusion  of  routines  to  compute  photodissociation  rates  and  rate  coefficients  with 
data  supplied  for  the  CFC-1 1  and  CFC-12  molecules,  and 

6.  Inclusion  of  user  control  of  solar  variability. 

These  changes  may  all  be  accessed  by  a  user  via  simple  flags  in  the  first  input  card.  The  model 
may  also  be  run  with  all  of  the  new  capabilities  turned  off,  providing  operation  identical  to  that  of 
either  MODTRAN  or  LOWTRAN7.  The  following  sub-sections  describe  each  AURIC 
modification. 


2.1  Schumann-Range  Absorption  Model 

The  Schumann-Runge  (SR)  bands  of  oxygen  represent  the  major  source  of  opacity  in  the 
important  UV  region  from  175  to  200  nm,^  and  are  an  important  source  of  O2  photodissociation.^ 
Recent  work  by  Minschwaner,  et  al,^’^  has  provided  an  accurate  SR  absorption  model  at  0.5  cm'^ 
resolution,  .4  polynomial  scheme  was  developed  by  Min.schwaner  and  the  Geophysics  Directorate 
to  approximate  the  temperature-dependence  of  the  spectral  cross-sections,  with  maximum  errors  of 
15%.^  The  algorithm  expresses  the  cross-section  as  a  quartic  polynomial  in  temperature,  with  only 
the  even  terms  non-zero.  The  spectral  dependence  resides  entirely  in  the  polynomial  coefficients. 
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This  formula  allows  the  optical  depth  to  be  factored  into  spectrally-dependent  and 
spectrally-independent  factors.  The  spectral  portion  consists  of  the  polynomial  coefficients, 
tabulated  every  0.5  cm’l  from  49,000  to  57,000  cm'^  For  accuracy,  separate  coefficients  were 
developed  for  each  of  three  temperature  regions,  from  130  to  190  K,  190  to  280  K,  and  280  to 
500  K.  The  spectrally-independent  portion  consists  of  only  nine  values,  namely  the  atmospheric 
path  column  density,  weighted  by  the  temperature  to  the  fourth,  second,  and  zeroth  powers, 
integrated  separately  in  each  of  the  three  temperature  ranges.  Such  temperature-weighted  densities 
were  previously  incorporated  into  MODTRAN  for  the  O3  Hartley-Huggins  band. 

MODTRAN  is  designed  to  step  spectrally  in  1  cm-I  increments,  with  the  output  optionally 
filtered  to  lower  resolution.  The  molecular  line  absorption  coefficients  are  stored  at  the  necessary 
1  cm-l  resolution,  while  slowly-varying  absorbers  (e.g.,  molecular  continua,  aerosol  extinctions) 
are  accessed  in  5  cm-1  steps  (usually  with  20  cm-l  resolution),  and  are  interpolated  to  the  1  cm-* 
wavenumber  grid. 

For  the  Schumann-Runge  model  of  Minschwaner,  et  al,  a  step  of  0.5  cm-*  is  desired. 
However,  major  changes  would  be  needed  to  modify  the  many  MODTRAN  routines  from  integer 
1  cm-*  steps  to  smaller  floating-point  steps.  G.  Anderson  and  L.A.  Hall  of  the  Phillips  Laboratory 
investigated  several  methods  of  reducing  the  Schumann-Range  resolution  to  1  cm-*,  and  concluded 
that  sufficient  accuracy  was  achieved  by  simply  sub-sampling  the  Minschwaner  coefficients,  taking 
every  other  one.^  A  routine,  "02SRBD",  was  written  to  read  in  the  Schumann-Range  spectral 
coefficients  within  the  desired  spectral  range,  and  to  return  the  values  for  any  wavenumber. 

The  non-spectral  portion  of  the  model  requires  the  path  column  density  of  O2,  scaled  by  the 
(normalized)  temperature  to  the  fourth,  second,  and  zeroth  powers,  and  stored  separately  for 
temperatures  in  each  of  three  temperature  ranges.  The  number  of  "species"  was  therefore  increased 
by  9,  and  subroutine  GEO  was  modified  to  compute  and  store  the  scaled  column  amounts. 

Although  the  SR  line  absorption  extends  up  to  57,000  cm-*,  the  line  wings  form  a  quasi¬ 
continuum  outside  the  band.  The  Minschwaner  model-^  provides  this  contribution,  which 
Minschwaner  has  fit  to  a  quadratic  polynomial  in  wavenumber,  from  57,000  to  66,000  cm-*. 

The  resolution  of  the  model  is  shown  in  the  spectral  U-ansmittance  plot  of  Figure  1 . 

2.2  Layering  Increase 

The  need  to  add  higher  altitudes  to  the  atmospheric  profiles  caused  the  need  for  an  increase 
in  the  maximum  number  of  layers  in  MODTRAN  from  the  original  34.  It  was  also  known  that  this 
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Figure  1 .  Sample  Transmittance  Calculations  In  the  Schumann-Runge 
Band  for  a  Vertical  Path  from  40  km  to  41  km 


modification  would  provide  benefits  in  other  ways,  and  at  all  wavelengths,  by  allowing  users  to 
utilize  more  than  34  layers  for  other  reasons.  For  example,  long  near-tangential  paths  are  less 
prone  to  errors  in  refraction  and  radiation  transfer  when  the  profile  is  finely  layered.  Having  a 
larger  dimension  also  releases  the  user  from  having  to  tailor  his  profile,  with  many  layers  at  some 
altitudes  and  very  few  layers  elsewhere.  Extra  layers  are  also  sometimes  desired  within  and  near 
inserted  cloud  decks. 

There  are,  actually,  several  different  types  of  layering  within  MODTRAN,  although  this  is 
not  obvious  to  the  user.  For  example,  there  are  "data  layers",  at  which  the  profile  is  input  or 
internally  stored,  and  "calculation  layers"  (with  parameters  interpolated  from  the  data  layers)  on 
which  the  computation  is  performed.  This  concept  is  a  good  one,  as  it  allows  for  calculations  to  be 
made  with  many  layers,  where  needed  for  accuracy,  while  keeping  the  input  layers  fixed. 

The  layering  schemes  currently  implemented  within  MODTRAN2  are: 

1 .  The  hard-coded  data  layering  for  Models  1  through  6,  set  in  block  data  MLATBM 
with  variable  ALT,  and  passed  in  common  block  /MLATM/. 
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2.  The  hard-coded  calculation  layering  for  Models  1  through  6,  set  in  subroutine 
FLA  YZ  with  variables  ZAER,  ZNEW,  and  ZNEWV,  and  used  locally  to  set  the 
array  ZMDL  for  use  elsewhere  in  the  program. 

3 .  The  user  input  data  layers  for  Model  7,  which  become  the  calculation  layers  in  that 
case. 

4.  Other  specific  data  layering  exists  for  clouds,  rain  rate,  aerosol  profile,  and  the 
Army  Vertical  Sumcture  Algorithm,  VSA.  These  are  independent  from  the  above 
layer  types,  and  are  inserted  into  the  calculation  layering  as  appropriate  by  routine 
FLAYZ. 

At  present,  no  independent  user  input  of  the  calculation  layering  is  allowed.  This  prevents  the  user 
from  independently  inputing  data  at  given  altitudes,  while  requesting  the  program  to  do  the 
calculations  at  different  (e.g.  more)  layers.  By  inputing  only  altitude  data  in  Model  7  with  the 
proper  flags  set,  a  user  is  able  to  run  with  Model  1  through  6  profiles,  calculated  at  the  input 
altitudes. 

Increasing  the  "34  layer"  dimension  was  a  large  task,  since  it  was  strongly  embedded 
("hardwired")  within  MODTRAN,  along  with  its  half,  double,  and  triple  (17,  68,  and  102).  Many 
variable  dimensions,  and  loop  limits,  contained  these  values  explicitly.  In  order  to  modify  these, 
and  to  make  future  changes  much  easier,  a  Fortran  PARAMETER  "LAYDIM"  was  defined  as  the 
maximum  number  of  layers  for  calculation,  and  was  increased  from  34  to  100.  The  multiples  of 
34  were  replaced  with  other  PARAMETERS,  caUed  LAYHAF,  LAYTWO,  and  LAYTHR.  All  of 
these  values  were  placed  in  a  single  file,  called  "AURIC.DEF",  which  is  automatically  inserted  in 
each  required  routine  via  the  Fortran  "INCLUDE"  command.  If  a  user  needs  to  increase 
LAYDIM,  he  need  only  modify  AURIC.DEF,  and  recompile  the  entire  program.  Although  the 
INCLUDE  command  is  not  ANSI  standard,  it  is  available  on  nearly  all  modem  compilers.  On 
other  compilers,  an  editor  must  be  used  to  explicitly  insert  the  file  in  place  of  the  INCLUDE 
statement. 

The  Air  Force  desired  to  allow  a  user  to  run  AURIC  in  a  mode  with  all  of  the  new  options 
turned  off,  including  the  layering.  For  this  reason,  the  original  calculation  layer  altitudes  were  kept 
in  place  (array  ZAER),  and  new  ones  were  added  in  a  separate  variable  (ZAUR),  used  only  when 
the  AURIC  mode  is  on.  The  ZAER  altitudes  vary  in  1  km  steps  from  0  to  25  km,  in  5  km  steps  up 
through  50  km,  with  other  layers  at  70,  100,  and  1000  km  (34  altitudes).  The  new  ZAUR 
altitudes  vary  in  1  km  steps  from  0  to  26  km,  in  2  km  steps  up  to  70  km,  in  4  km  steps  up  to 
150  km,  in  10  km  steps  up  to  300  km,  with  other  layers  at  500, 750,  and  1000  km  (87  altitudes). 
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In  MODTRAN,  the  default  ZAER  altitudes  are  modified  if  the  user  selects  one  of  several 
different  aerosol  algorithms,  including  the  VS  A  and  cloud  models.  These  options  are  also 
supported  in  AURIC  mode,  with  the  appropriate  aerosol  model  altitudes  replacing  the  ones  in 
ZAUR. 

2.3  Altitude  Extension 

MODTRAN  contains  six  model  atmosphere  profiles  within  it,  contained  in  block  data 
MLATMB.  There  are  50  "data  layer"  altitudes  defined  (the  number  is  now  set  in  AURIC  via  the 
Fortran  parameter  MODDAT).  These  profiles  step  by  1  km  from  0  to  25  km,  by  2.5  km  up  to 
50  km,  and  by  5  km  up  to  120  km.  These  data  layers  are  still  accessed  when  AURIC  is  run  in 
MODTRAN  mode. 

Two  new  sets  of  data  layers  were  added  for  AURIC,  one  for  daytime  (block  data 
ALATMB)  and  one  for  nighttime  (block  data  NLATMB).  In  the  altitude  regime  from  0  to  50  km, 
the  parameters  were  taken  from  the  MODTRAN  data  set,  but  interpolated  onto  a  different  altitude 
grid.  From  50  to  300  km,  parameters  from  the  SHARC  model^  were  used  for  Model  6,  and  zeros 
elsewhere.  Above  300  km,  values  from  the  1976  U.S.  Standard  Atmosphere  were  used  for 
Model  6,  and  zeros  elsewhere.  A  total  of  104  layers  were  defined,  with  altitudes  from  0  to  25  km 
by  1  km,  up  to  50  km  by  2.5  km,  up  to  150  km  by  2  km,  up  to  300  km  by  10  km,  with  additional 
layers  at  500, 750,  and  1000  km.  Models  1  through  5  were  not  filled  out  above  50  km  because  of 
the  lack  of  appropriate  SHARC  and  1976  U.S.  Standard  data.  In  addition,  the  U.S.  Standard 
Aunosphere  lacks  data  above  300  km  for  all  components  except  N2,  O2,  and  the  total  air  density. 

2.4  Solar  Irradiance  Wavelength  Extension 

The  extraterrestrial  solar  irradiance  data  stored  within  the  model  were  extended  to  shorter 
wavelengths,  using  the  1987  SUSIM  data  taken  from  Spacelab2.^  The  MODTRAN  data  extend  up 
to  57,470  cm‘f  The  new  data  extend  from  25,003  to  82,987  cm’^  (120.5  to  399.95  nm)  at 
0.05  nm  spacing,  and  with  0. 15  nm  spectral  resolution.  The  new  data  are  used  above  25,000  cm-^ 
in  AURIC  mode.  In  MODTRAN  mode,  the  original  data  set  is  used. 

Since  the  SUSIM  data  are  tabulated  versus  wavelength  (nm),  routine  SUN  was  modified  to 
access  the  appropriate  data  set,  and  to  interpolate  in  wavelength  where  necessary. 
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2.5  Photodissociation  Computation 


The  photodissociation  of  many  molecules,  including  that  of  chloro-fluoro-carbons  (CFCs), 
depends  on  the  absorption  of  sunlight  by  oxygen  in  the  Schumann -Range  and  Herzberg  bands. 
Since  AURIC-M  contains  an  improved  Schumann-Range  model  it  is  a  natural  vehicle  for  the 
computation  of  photodissociation  ("PD")  rates  and  rate  coefficients.  AURIC-M  was  modified  to 
make  these  computations  for  the  CFC-1 1  (CFCI3)  and  CFC- 12  (CF2CI2)  molecules  as  a 
demonstration  of  the  model's  capabilities.  Other  molecules  may  be  easily  added  to  the  database. 

Changes  to  AURIC-M  were  made  in  three  areas: 

•  The  addition  of  CFC  concentration  profiles  to  the  internal  data  structures; 

•  The  generation  of  a  data  file  containing  CFC  photodissociation  cross-sections 
versus  temperature  and  wavelength;  and, 

•  The  creation  of  a  "PD  mode",  in  which  the  PD  rates  are  computed,  instead  of 
U'ansmittance  and  radiance. 

Each  of  these  topics  are  discussed  below. 

2.5.1  Addition  of  CFC  Concentration  Profiles 

Mixing  ratios  for  CFCs  had  already  been  incorporated  into  a  developmental  version  of 
MODTRAN,  the  model  on  which  AURIC-M  is  based.  Unfortunately,  the  original  AURIC 
development  changed  all  of  the  internal  molecular  storage,  so  that  it  has  duplicate  storage  for  both 
the  MODTRAN  and  AURIC-M  modes  (see  section  2.2).  The  MODTRAN  CFC  changes  thus 
needed  to  be  added  into  AURIC-M  twice,  in  parallel.  A  full  integration  was  done  so  that  the  CFC 
profiles  are  accessible  to  the  entire  model,  including  the  ability  to  use  both  a  hard-coded 
("MODEL=6")  profile  and  user-defined  ("MODEL=7")  inputs.  This  integration  involved  changing 
routines  DRIVER,  BMDATA,  DEFALT,  AERNSM,  STDMDL,  and  VSANSM.  Although  an 
easier  implementation  could  have  been  done,  by  using  only  a.  hard-coded  profile  accessible  only  by 
the  PD  routines,  the  current  implementation  has  several  advantages,  including: 

•  It  allows  for  other  types  of  CFC  computation  by  other  portions  of  the  AURIC-M 
code; 

•  It  retains  a  closer  relationship  between  AURIC-M  and  MODTRAN;  and, 

•  It  is  much  more  flexible  for  the  user. 

CFC  concentration  profiles  were  provided  by  the  Phillips  Laboratory,  in  the  Fortran  block 
data  format  used  for  other  molecules,  using  the  MODTRAN  altitude  grid  (see  Sections  2.2  and 
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2.3).  These  data  were  then  exponentially  interpolated  onto  the  larger  AURIC  altitude  grid.  The 
MODTRAN  data  were  provided  in  a  block  data  routine  XMLATM.  The  "X"  designation  is 
because  they  are  called  "cross-section"  molecules  in  the  supplied  routines  (molecules  with 
absorptions  given  as  cross-sections,  and  not  as  line  parameters).  We  have  kept  the  "X" 
designation  on  most  of  the  related  variables.  The  AURIC  interpolated  version  of  these  profiles  are 
defined  in  the  new  block  data  XALATM.  Although  storage  is  set  up  for  18  such  species,  profile 
data  is  supplied  for  the  following  eight:  CFC-1 1,  CFC-12,  CFC-13,  CFC-14,  CFC-22,  CFC-1 13, 
CFC-114,  and  CFC-1 15. 

The  standard  MODTRAN/AURIC  molecule  concentrations  pass  through  several  variables 
(AMOL,  WMOL,  etc.),  in  units  of  molecules/cm-^,  as  they  are  processed  by  routines  DRIVER, 
AERNSM,  DEFALT,  and  RDNSM.  Parallel  variables  were  defined  for  the  cross-section 
molecules  in  these  routines.  The  concentrations  are  then  converted  by  routine  STDMDL  to  units 
proportional  to  the  standard  spectroscopic  unit  of  "amagat",  and  are  stored  in  array  DENSTY 
versus  absorption  "band"  (Fortran  parameter  MBANDS)  instead  of  just  per  molecule.  For  the 
cross-section  species,  the  same  variable  DENSTY  is  used,  except  that  the  array  size  MBANDS 
was  increased  to  hold  the  extra  storage.  From  there  on,  no  code  changes  were  needed  to  store  the 
cross-section  molecules,  except  for  the  increase  in  the  value  of  MBANDS.  After  the  calls  to 
routine  GEO  the  amounts  of  all  molecules,  integrated  over  the  various  atmospheric  paths,  are 
available  in  arrays  WPATH  and  W,  in  units  of  amagat-cm. 

2.5.2  Photodissociation  Cross-Section  Database 

Molecules  CFG- 1 1  and  CFC-12  were  selected  for  initial  input  into  the  cross-section 
database  because  of  their  high  level  of  current  research  interest,  and  the  ready  availability  of  both 
concentration  profile  and  cross-section  data  for  them.  The  cross-sections  were  taken  from 
DeMore,  et  al.,2  where  they  are  tabulated  versus  wavelength  for  a  temperature  of  298  K. 
Expressions  are  also  given  for  the  variation  with  temperature.  Figure  2  shows  the  temperature¬ 
scaling  parameter  used. 

The  data  were  put  into  an  ASCII  datafile  which  is  designed  to  have  more  species  added 
easily.  The  analytic  temperature  dependence  was  replaced  with  tabulated  values  computed  from  the 
given  expression  in  order  to  keep  the  database  as  generic  as  possible  (i.e.,  in  order  to  not  assume  a 
specific  temperature-dependence).  A  new  routine  called  PDCRSS  was  written  to  read  that  datafile, 
and  to  interpolate  the  data  on  temperature  and  wavelength. 
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F11  Temperature  Scaling 


t(x.y)  — 


Figure  2.  Scaling  Factor  for  the  CFCl  1  Absorption  Cross-Section  Versus 
Wavelength  and  Temperature 

2.5.3  Photodissociation  Rate  Computations 

MODTRAN  and  AURIC-M  have  always  had  the  ability  to  perform  several  modes  of 
computation,  controlled  through  an  input  variable  "lEMSCT".  The  modes  vary  from  transmittance 
(IEMSCT=0),  through  thermally-emitted  radiance  (IEMSCT=1),  to  solar-scattered  radiance 
(IEMSCT=2).  Computations  of  incident  solar  irradiance  are  computed  with  IEMSCT=3.  For  this 
effort  a  new  "PD"  mode  was  defined,  controlled  by  setting  lEMSCT  to  4. 

A  photodissociation  rate  "L",  in  [l/s/cm^],  is  computed  as: 


L  =  JN 


(1) 


where  N  is  the  molecular  number  density  [  l/cm^]  and  J  is  the  photodissociation 
rate  coefficient.  The  coefficient!  is  computed  from  the  molecule’s 
cross-section  and  the  incident  irradiance: 


J  = 


dv 


aE 

(4 


(2) 
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where 


V  is  wavenumber  (cm'^), 

<7  is  the  wavenumber-dependent  PD  cross-section  (cm^/photon), 

(hv)  is  the  energy  per  photon  [J/photon],  and 
Ey  is  the  incident  spectral  irradiance  [W/cm^/cm’^]. 

The  computations  of  N  and  o  were  described  previously.  The  other  required  computation  is 
obviously  the  solar  irradiance  incident  at  an  altitude.  Computations  at  multiple  altitudes  are  also 
desired. 

The  solar-scatter  mode  (IEMSCT=2)  of  AURIC-M  already  computes  the  solar  irradiance 
incident  on  each  altitude  layer  of  a  sensor  line-of-sight  (LOS).  By  setting  up  a  fictitious  vertical 
LOS  between  two  altitudes  of  interest,  the  PD  computation  can  piggy-back  on  the  existing 
algorithm  to  compute  the  desired  irradiance,  at  each  layer  altitude  along  the  specified  LOS  path.  A 
special  PD  routine  is  then  called  to  compute  the  PD  rate  for  each  of  these  altitudes. 

The  first  step  in  the  above  procedure  is  to  initialize  the  software.  When  the  user  calls  for 
the  PD  mode  by  setting  IEMSCT=4,  several  variables  are  set: 

A  new  logical  variable  LPHDIS  is  set  to  TRUE,  for  later  testing  on  whether  or  not 
the  PD  mode  is  set; 

lEMSCT  is  reset  to  2,  to  allow  the  existing  algorithms  to  compute  the  required  solar 
irradiances; 

IMULT  is  set  to  0,  to  force  the  neglect  of  multiple  scattered  solar  calculations; 

Variable  LAURIC  is  set  to  TRUE,  to  force  AURIC  mode  on; 

Variable  ITYPE  is  set  to  2,  to  force  the  LOS  path  geometry  to  be  a  slant  path 
between  two  altitudes; 

Variable  THETA  is  set  to  zero,  to  force  the  LOS  path  to  be  a  vertical  one;  and. 

Variables  H 1  and  H2  are  forced  to  be  in  increasing  order,  so  that  the  LOS  path  start 
and  end  altitudes  are  consistent  with  a  value  of  THETA=0. 

Note  that  any  values  of  IMULT,  LAURIC,  ITYPE,  and  THETA  input  by  the  user  are  over-ridden 
by  the  above.  Multiple  solar  scatter  is  neglected  because  of  the  algorithmic  complexity  of  including 
it.  In  addition,  computations  by  Minschwanerl®  have  shown  that  the  effect  is  small  (light  scattered 
in  is  mostly  offset  by  light  scattered  out). 
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The  computation  of  incident  sunlight  is  performed  in  two  parts.  Existing  routine  LOOP 
(called  for  each  wavenumber  by  routine  TRANS)  computes  the  transmittance  from  the  sun  to  each 
altitude  along  the  LOS  path.  In  the  normal  solar-scatter  mode  this  transmittance  is  combined  with 
that  along  the  LOS  path  to  the  observer.  LOOP  was  therefore  modified  so  that,  in  PD  mode,  only 
the  transmittance  to  the  LOS  is  included,  and  not  the  transmittance  along  it.  The  set  of  solar 
transmittances  to  each  altitude  are  then  passed  back  to  routine  TRANS,  which  calls  a  new  routine 
PHTDIS  to  perform  the  actual  rate  computation. 

PHTDIS  has  inputs  of  the  molecular  PD  cross-sections  (from  PDCRSS),  the  solar 
transmittances  (from  LOOP),  the  current  wavenumber  (from  TRANS),  and  the  molecular 
concentrations  per  altitude  (from  GEO).  It  calls  existing  routine  SOURCE  to  obtain  the  extra¬ 
terrestrial  solar  irradiance,  multiplying  that  by  the  solar  transmittances.  It  then  computes  the  PD 
rates  and  rate  coefficients  from  Eqs.  (1)  and  (2),  replacing  the  integral  with  a  sum  over 
wavenumber.  The  results  are  returned  to  TRANS  for  output,  versus  altitude  and  molecular 
species. 

A  sample  computation  is  shown  in  Figure  3,  for  the  standard  atmosphere  and  the  sun  at  60° 
zenith.  Figure  3a  plots  the  rate  coefficients  for  CFCl  1  and  CFC12  versus  altitude,  and  Figure  3b 
plots  the  rates. 

2.6  Solar  Variability 


The  solar  optical  flux  is  known  to  vary  over  time,  with  a  dominant  cycle' '  of  1 1  years  and 
a  smaller  cycle  of  27  days.  Brasseur"  reviewed  the  available  data  and  generated  a  table  of  peak- 
to-peak  variations  (in  percent)  versus  wavelength,  from  1 16  to  303  nm  for  the  1 1-year  cycle.  He 
also  uses  one-half  those  variations  for  the  27-day  cycle.  His  1 1-year  values  are  shown  in 
Figure  4. 

For  AURIC-M  the  level  of  solar  activity  is  expressed  by  the  user  in  terms  of  two  new  input 
variables  al  1  and  a27.  Each  of  these  denotes  the  relative  level  of  the  corresponding  cycle,  with  a 
value  of  zero  meaning  the  cycle  minimum  and  one  meaning  the  maximum.  Given  a  minimum  solar 
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Figure  3.  Sample  Computation  of  Photodissociation  (a)  Rate  Coefficients 
and  (b)  Rates  for  CFCl  1  and  CFC12 
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Figure  4.  Solar  Relative  Variability  Versus  Wavenumber 


irradiance  at  a  wavenumber,  with  a  corresponding  maximum  level  of  E^ax’  we  can  write  (for 
each  cycle)  an  arbitrary  irradiance  E  in  terms  of  "a"  by: 


a 


E-E 


min 


E  -E 
max  min 


(3) 


Brasseur's  tabulated  peak-to-peak  relative  variations  can  then  be  expressed  as: 


E  -  E 

max  min 

E  . 

mm 


(4) 


Combining  these  two  equations  provides: 


E  —  E|uj|]  ( 1  +  aV)  .  (5) 

If  we  now  assume  that  the  two  cycles  are  independent,  and  additive,  we  can  derive  an  expression 
for  the  irradiance  for  given  al  1  and  a27  values  as: 

E(al  l,a27)  =  E^in  (I  +  al  IVl  1  +  a27V27)  ,  (6) 
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where  VI 1  and  V27  are  the  variation  levels  from  Brasseur’s  tabulation.  Note  that  E,  Emin’ 
and  V27  all  vary  with  wavelength. 

The  AURIC-M  input  was  modified  to  read  in  the  values  of  al  1  and  a27  when  the 
photodissociation  mode  is  turned  on.  These  are  used,  along  with  stored  spectral  values  for  the 
solar  variation  and  the  "baseline"  solar  irradiance,  to  generate  the  desired  irradiance  from  Eq.  (6). 

3.0  MODEL  USAGE  CHANGES 

Most  of  the  changes  from  MODTRAN  to  AURIC-M  were  designed  to  not  require  any 
additional  user  input.  A  single  input  flag  (described  below)  turns  the  AURIC  mode  on  or  off. 
Several  new  data  files  are  required,  but  the  user  need  not  modify  them. 

Section  3. 1  describes  the  AURIC-M  inputs,  and  comments  on  the  user  restrictions  in 
altitude  profiles  and  spectral  parameters.  Section  3.2  describes  the  new  data  files. 

3.1  Changes  To  The  MQDTRAN2  User  Input 

MODTRAN  was  designed  to  operate  in  two  primary  modes,  either  using  the  MODTRAN 
capabilities  (higher  spectral  resolution  and  temperature-dependent  band  models),  or  defaulting  to 
the  operation  of  LOWTRAN7,  the  model  on  which  MODTRAN  was  built.  The  run  mode  is 
controlled  with  a  single  input  "flag"  variable.  While  the  MODTRAN  mode  is  more  accurate,  it 
runs  more  slowly,  and  does  not  extend  spectrally  beyond  22,681  cm'^  (440  nm). 

AURIC  has  similarly  been  designed  with  the  user  option  for  either  LOWTRAN7, 
MODTRAN,  or  AURIC  mode.  The  primary  user  input  change  from  MODTRAN  is  the 
substitution  of  a  three-space  character  variable,  "TYPROG",  for  the  previous  MODTRAN  logical 
variable.  MODTRAN  is  turned  on,  rather  than  LOWTRAN  7,  if  a  "T"  character  is  in  the  first 
space  (for  compatibility  with  MODTRAN).  MODTRAN  is  allso  turned  on  if  an  "M"  is  in  any  of 
the  three  spaces.  AURIC  is  turned  on,  independent  of  MODTRAN/LOWTRAN,  if  an  "A"  is  in 
any  of  the  three  spaces.  Since  the  upper  atmosphere  profile  is  day-  and  night-dependent,  an  "N" 
character  in  any  of  the  three  spaces  will  be  used  with  AURIC  to  specify  night.  The  characters  "L" 
and  "D"  may  also  be  used  to  explicitly  indicate  LOWTRAN  and  "day"  in  the  absence  of  "M"  and 
"N".  Note,  however,  that  all  of  the  above  character  flags  must  be  upper  case  letters.  The 
controlling  Fortran  statements  in  subroutine  DRIVER  are: 
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LOGICAL  MODTRN,AURICL, NIGHT 
CHARACTER* 3  TYPROG 
COMMON  /PRGTYP/  AURICL, NIGHT 
MODTRN=. FALSE. 

AURICL=. FALSE. 

NIGHT  =. FALSE. 

READdRD,  •  (A3,I2,12I5,F8.3,F7.2)  '  )  TYPROG,  MODEL,  ITYPE,  lEMSCT, 
&  IMULT , Ml , M2 , M3 , M4 , M5 , M6 , MDEF , IM , NOPRT , TBOUND , SALE 
IF (TYPROG (1:1) . EQ . ' T ' )  MODTRN= . TRUE . 

IF (TYPROG ( 1 : 1 ) . EQ . ' M ■ )  MODTRN= . TRUE . 

IF (TYPROG(2 :2) .EQ. ’M' )  MODTRN= . TRUE . 

IF (TYPROG(3 :3) .EQ. ’M' )  MODTRN= . TRUE . 

IF (TYPROG (1: 1) .EQ. 'A‘ )  AURICL= . TRUE . 

IF (TYPROG (2: 2) .EQ. ’A' )  AURICL= . TRUE . 

IF(TYPROG(3:3) .EQ. 'A' )  AURICL= . TRUE . 

IF (TYPROG (1: 1) .EQ. 'N‘ )  NIGHT  =.TRUE. 

IF(TYPROG(2:2) .EQ. ’N' )  NIGHT  =.TRUE. 

IF  (TYPROGO  :3)  .EQ.  'N' )  NIGHT  =.TRUE. 

3.1.1  Altitude  and  Wavelength  Options 


AURIC  is  delivered  with  the  maximum  number  of  calculational  layers  set  to  a  value  of 
LAYDIM  =  100.  This  may  be  changed  by  editing  the  file  AURIC.DEF,  which  gets  included  in  the 
source  code  via  the  compiler.  The  other  calculational  parameters  (e.g.,  LAYTWO)  are  defined  via 
scaling  LAYDIM,  so  the  user  need  not  modify  them.  The  code  must  then  be  re-compiled  for  the 
re-dimensioning  to  take  place. 


The  user  may  also  wish  to  modify  the  default  calculational  layers  (used  except  when 
Model=7),  which  are  defined  in  array  ZAUR  within  subroutine  FLAYZ.  These  87  layers  are  those 
used  with  atmospheric  Models  1  through  6.  If  they  are  not  modified,  the  user  can  only  calculate 
with  more  layers  by  the  use  of  Model  7.  Runs  with  the  six  standard  models  will  use  the  87 
calculation  layers  in  ZAUR,  plus  any  others  called  for  by  the  user's  selection  of  special  aerosol 
options,  including  the  Vertical  Structure  Algorithm  (VS A)  and  cloud  layers.  These  options  can 
currently  result  in  the  use  of  up  to  96  calculational  layers. 


Since  the  user-defined  Model  7  allows  all  of  the  parameters  (temperature,  pressure, 
concentrations)  to  be  interpolated  from  another  model  (e.g.,  from  Model  6)  via  the  variables  Ml 
through  M6  plus  MDEF,  the  user  may  define  a  Model  7  input  which  specifies  only  the  desired 
altitudes,  and  lets  AURIC  insert  the  proper  parameters  from  a  standard  model  (1  to  6)  at  each 
altitude.  Note  that  when  Model=7,  the  calculational  layers  are  the  same  as  the  input  data  layers.  If 
the  default  calculational  layers  in  ZAUR  are  not  changed,  the  user  will  always  get  the  same  ones 
when  Models  1  through  6  are  used. 
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The  existing  MODTRAN  code  allows  a  user  to  input  his  own  profile,  through  the 
"MODEL=7"  option.  In  this  situation  input  cards  "2C"  and  "2C 1  must  be  supplied.  If  variable 
IRD1=1,  then  card  "2C2"  is  also  read,  containing  concentrations  of  eight  molecules,  including 
N2O,  CO,  and  CH4.  A  card  "2C2B"  is  then  read  with  the  concentration  of  HNO3.  Eight  new 
molecules  were  added  to  AURIC-M  for  the  photodissociation  option.  These  "cross-section" 
molecules  are:  CFC-1 1,  CFC-12,  CFC13,  CFC-14,  CFC-22,  CFC-1 13,  CFC-1 14,  and  CFC- 
115.  When  MODEL=7  and  IRD1=1,  concentrations  for  these  eight  molecules  will  be  read  from  a 
new  card  "2C2C",  directly  following  card  2C2B  at  each  altitude.  The  user  has  the  same  options 
for  input  units  with  these  as  he  does  for  the  other  molecules,  which  are  controlled  by  the  14 
"JCHAR"  variables  read  in  on  card  2C 1 .  The  CFC  units  are  set  by  a  new  variable  JCHARX,  read 
in  on  card  2C1  immediately  following  the  other  14  JCHARs.  Note  that  the  single  JCHARX  value 
controls  the  units  for  all  of  the  CFCs. 


The  new  MODEL=7  format  is  thus  as  follows: 


CARD  2C  ML,  IRD 1 ,  IRD2,  TITLE 

FORMAT  (315,  18A4) 

CARD  2C1  ZMOL,P,T,WMOL(l),  WMOL(2),  WMOL(3),  (JCHAR(J),  J  >  1,  14),  JCHARX 
FORMAT  (FIO.3;  5E10.3,  15A1) 


CARD  2C2  (WMOL(J),  J  =  4, 1 1 ) 
FORMAT  (8E10.3) 

CARD  2C2B  WMOL(  1 2) 

FORMAT  (8E  10.3) 

CARD  2C2C  WMOLX(J),  J  =  1 ,8) 
FORMAT  (8E10.3) 


[ONLY  IF  IRD  1  =  1] 
[ONLY  IF  IRD  1  =  1] 
[ONLY  IF  IRD  1  =  1] 


Cards  2C3  and  later  are  the  same  as  for  MODTRAN. 


MODTRAN  is  restricted  to  wavenumbers  below  50,000  cm'^  (wavelengths  above 
200  nm),  and  only  the  LOWTRAN7  mode  is  available  from  22,68 1  to  50,000  cm'^  (200  to 
440  nm).  In  the  LOWTRAN7  mode,  the  computations  are  made  at  20  cm'*  spectral  resolution,  in 
step  sizes  which  are  a  multiple  of  5  cm'*.  In  the  MODTRAN  mode  the  computations  are  made 
internally  in  1  cm'*  steps,  and  are  spectrally  degraded  to  a  user-specified  resolution  and  step  size. 
In  the  AURIC  mode,  the  spectral  step  size  is  dependent  on  the  spectral  range,  and  on  whether  or 
not  the  MODTRAN  flag  was  set.  When  the  entire  spectral  range  is  in  the  UV  region  (above 
25,000  cm'*;  below  400  nm)  the  internal  step  size  is  set  to  1  cm**,  and  the  output  is  degraded  as 
requested  by  the  user.  When  part  of  the  spectral  range  is  in  the  infrared  or  visible,  the  spectral 
resolution  is  set  by  the  MODTRAN  flag  (i.e.,  1  cm  *  steps  if  MODTRAN  is  on,  or  a  multiple  of 
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5  cm'^  if  not.)  Table  1  summarizes  the  utilization  of  various  options  versus  the  run  mode.  It 
shows  some  of  the  interplay  between  the  MODTRAN  and  AURIC  flags,  coupled  with  the 
wavenumber  restrictions  on  them. 


Table  1 .  Status  of  AURIC  Layering  and  Spectral  Options  Versus  Run  Mode 


Mode  Name 

MODTRAN 
Flag  On? 

AURIC 
Flag  On? 

Starting 

Wavenumber 

(cm'I) 

Ending 

Wavenumber* 
(cm- 1) 

Internal 
Step  Size 
(cm-l) 

Maximum 

#of 

Layers 

LOWTRAN? 

N 

N 

COi  <  CD2 

any 

N*5b 

34 

MODTRAN 

Y 

N 

CO]  <C02 

0)2  <  22,681 

1 

34 

LOWTRAN7C 

Y 

N 

0)1  <0)2 

0)2  >  22,681 

N*5c 

34 

AURIC/IR-Vis 

N 

Y 

0)]  <  25,000 

CO]  <  CO2 

N*5 

100 

AURIC/UV 

N 

Y 

0)]  >  25,000 

CO]  <  CO2 

Id 

100 

AUR/MOD/IR-Vis 

Y 

Y 

0)1  <  W2 

C02<  22,681 

1 

100 

AUR/LOW 

Y 

Y 

0)i  <  25,000 

CO2  >  22.68 1 

N*5 

100 

AURIC/UV 

Y 

Y 

®l  ^25,(X)0 

CO]  <  CO2 

Id 

100 

a  All  wavenumbers  must  be  less  than  83,000 
b  Means  any  positive  integer  multiple  of  5 

c  Beyond  band  model  data;  MODTRAN  mode  is  turned  off  internally 
d  For  Schumann-Runge  Model _ 

3.1.2  Photodissociation  Model 


In  order  to  perform  photodissociation  rate  computations  with  AURIC-M  the  variable  lEMSCT 
must  be  set  to  4.  This  will  automatically  put  the  run  into  AURIC  mode,  regardless  of  the  input  value 
of  TYPROG  (described  above).  Other  variables  automatically  set  are; 


•  IMULT  =  0; 

•  ITYPE  =  2; 

•  THETA  =  0.0; 
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the  input  values  for  these  will  be  over-ridden  by  the  above  values.  The  other  inputs  required  of  the 
user  are: 

•  The  atmospheric  profile; 

•  The  start  and  end  altitudes,  HI  and  H2; 

•  The  solar  position;  and, 

•  The  spectral  range. 

The  atmospheric  profile  inputs  are  the  same  as  for  MODTRAN,  with  the  exception  of  the  MODEL=7 
case,  as  described  in  Section  3.1.1,  when  the  additional  JCHARX  and  WMOLX  inputs  are  read  on 
cards  2C1  and  2C2C.  Note  that  the  current  PD  cross-section  database  only  includes  molecules  CFC- 
1 1  and  CFC-12  (the  first  two  columns  of  Card  2C2C)  -  concentrations  for  the  other  six  cross-section 
molecules  thus  do  not  matter. 

The  PD  computation  forces  the  pseudo  "sensor  look  geometry"  into  a  vertical  path 
(THETA=0.0)  defined  by  the  two  altitudes  HI  and  H2.  AURIC-M  automatically  forces  HI  <  H2,  so 
the  user  doesn't  need  to.  Variables  RANGE  and  BETA  will  not  be  used. 

The  solar  position  is  defined  as  with  MODTRAN,  including  the  various  "IP ARM"  input 
options.  Note  that  the  variable  lEMSCT  is  internally  reset  to  a  value  of  2,  so  that  the  solar  position  is 
required.  For  the  PD  case  the  user  must  also  input  the  values  of  the  solar  activity  levels  al  1  and  a27 
(see  section  2.6).  These  are  entered  on  a  new  Card  3A3,  directly  following  Card  3A2.  The  format  is: 

CARD  3 A3  A11,A27 
FORMAT  (2F5.3) 

Note  that  this  card  is  needed  if  and  only  if  the  user  entered  IEMSCT=4. 

3.2  Additional  Data  Files 

Three  types  of  fixed  data  files  have  been  added  to  the  model.  Their  formats  are  documented 

below. 

3.2.1  Schumann-Runge  Coefficient  Files 

Subroutine  02SRBD  reads  and  returns  the  polynomial  coefficients  for  the  O2  Schumann- 
Runge  band,  called  A,  B,  and  C.  Absorption  cross  sections  are  computed  from  these,  and  the 
temperature,  by: 
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Cross  section  (cm^)  =  ( A*X'^  +  B*X  +  C  )*  lO'^O , 


(7) 


where. 


X  =  [(100-T)/10]2, 

and  T  is  the  temperature  in  Kelvins.  Each  coefficient  is  tabulated  over  the  wavenumber  range  from 
49,001  to  57,000  cm'^  at  intervals  of  1  cm'^  and  over  the  three  temperature  ranges:  130K  to  190K, 
190K  to  280K,  and  280K  to  500K.  The  three  files  130190.COF,  190280.COF,  and  280500.COF  are 
read  the  first  time  through;  those  coefficients  within  the  user's  spectral  range  are  stored. 


Each  file  begins  with  three  header  lines,  which  are  not  used  by  AURIC.  A  list  of 
wavenumbers,  and  coefficients  for  each,  is  read  next.  The  file  format  is: 

Line  Type  #  Format  Variables  Type  Definition 

LINE  TYPE  1  A80 

LINE  CHAR*80  -  not  used  - 

-  Repeat  line  type  1  for  a  total  of  three  lines  - 
LINE  TYPE  2  F8.1,3E1 1.3,F6.2,F8.1 


WN 

REAL 

Wavenumber  (cm-l) 

TCA 

REAL 

Coefficient  A 

TCB 

REAL 

Coefficient  B 

TCC 

REAL 

Coefficient  C 

DUMPl 

REAL 

-  not  used  - 

DUMP2 

REAL 

-  not  used  - 

-  Repeat  line  type  2  for  a  total  of  8000  lines  - 
3.2.2  SUSIM  Solar  Irradiance  Data 

Subroutine  RDSUSM  reads  the  data  file  'IRRAD_HI.INP',  containing  the  high-resolution 
SUSIM  solar  irradiance  data.  This  data  set  covers  the  spectral  range  from  120.5  to  399.95  nm  in  steps 
of  0.05  nm  (5590  values).  The  original  data  file  had  the  irradiances  stored  10  per  line,  with  a 
wavelength  beginning  each  line.  Each  of  the  original  lines  has  been  divided  into  two,  so  that  each  line 
is  now  less  than  80  characters  wide,  for  easy  portability.  The  first  line  has  the  wavelength  and  six 
irradiance  values,  followed  by  the  other  four  values  on  the  next  line.  The  format  is: 
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Line  Tvpe  # 

Format 

Variables 

Type 

Definition 

LINE  TYPE  1 

1X,F6.2,6E12.4 

LAMBDA 

REAL 

Wavelength  (nm)  [not  used] 

IRRAD(J) 

REAL 

Irradiance  at  LAMBDA 

IRRAD(J+1) 

• 

• 

REAL 

Irradiance  at  LAMBDA  +  0.05nm 

• 

IRRAD(J-i-5) 

REAL 

Irradiance  at  LAMBDA  +  0.25nm 

LINE  TYPE  2 

E1L4,3E12.4 

IRRAD(J+6) 

• 

• 

REAL 

Irradiance  at  LAMBDA  +  0.30nm 

• 

IRRAD(J+9) 

REAL 

Irradiance  at  LAMBDA  +  0.45nm 

-  Repeat  the  above  pair  of  lines  a  total  of  559  times  - 

This  format  produces  a  file  with  1,118  lines  (559  line  pairs). 


3.2.3  Photodissociation  Cross-Section  Data 

Subroutine  PDCRSS  reads  the  file  "PDCROSS.INP",  containing  the  PD  absorption  cross- 
sections  versus  wavelength  and  temperature.  The  current  data  set  contains  values  for  the  molecules 
CFCl  1  and  CFC12,  covering  the  spectral  range  of  170  -  260nm,  and  temperatures  from  200  -  320K. 

The  file  format  is  straightforward.  The  first  line  gives  the  number  of  molecules  contained  in 
the  file.  All  of  the  data  for  the  first  molecule  is  given  next,  followed  by  all  of  the  data  for  the  second 
molecule,  etc.  For  any  given  molecule  the  number  of  temperature  points  is  listed,  followed  by  the 
temperature  values.  The  number  of  wavelength  points  is  then  listed,  followed  by  a  set  of  lines,  one 
line  per  wavelength.  Each  such  line  contains  the  wavelength  value,  followed  by  the  cross-sections  at 
that  wavelength,  for  each  temperature.  The  next  molecule  (if  any)  is  then  read. 

The  file  format  also  allows  for  cross-sections  which  are  independent  of  temperature.  In  this 
case  the  number  of  temperatures  is  input  as  one  or  less,  and  the  line  of  temperatures  is  skipped.  Note, 
therefore,  that  in  order  to  interpolate  in  temperature  at  least  two  temperature  values  must  be  given  in  the 
file. 
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The  format  is: 


Line  Type  #  Format  Variables  Type  Definition 

LINE  TYPE  1  15  NENTRY  INTEGER  Number  of  molecules  in  the  file 

-  Repeat  Line  Types  2  through  6  for  each  of  NENTRY  molecules 

LINE  TYPE  2  A  ENNAME  CHAR*8  Molecule  name 

LINE  TYPE  3  15  NTEMPS  INTEGER  Number  of  temperatures  (max  of  10) 

-  If  NTEMPS  >=  2  then  read  Line  Type  4 

LINE  TYPE  4  12E10.3  TEMPS(I)  REAL  Array  of  NTEMPS  temperatures  (K) 

LINE  TYPE  5  15  NWAVLS  INTEGER  Number  of  wavelengths  (max  of  1 00) 

-  Repeat  Line  Type  6  NWAVLS  times 

LINE  TYPE  6  12E10.3  WAVES  REAL  Wavelength  (nm) 

PDCR(I)  REAL  Array  of  NTEMPS  cross-sections  at  the 

current  wavelength 

-  Next  molecule  entry  (if  any) 
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